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Abstract 

As  the  first  gravity  wave  (GW)  climatology  study  using  nadir-viewing  infrared  sounders, 
50  Atmospheric  Infrared  Sounder  (AIRS)  radiance  channels  are  selected  to  estimate 
GW  variances  at  pressure  levels  between  2-100  hPa.  The  GW  variance  for  each  scan 
5  in  the  cross-track  direction  is  derived  from  radiance  perturbations  in  the  scan,  indepen¬ 
dently  of  adjacent  scans  along  the  orbit.  Since  the  scanning  swaths  are  perpendicular 
to  the  satellite  orbits,  which  are  inclined  meridionally  at  most  latitudes,  the  zonal  compo¬ 
nent  of  GW  propagation  can  be  inferred  by  differencing  the  variances  derived  between 
the  westmost  and  the  eastmost  viewing  angles. 
io  Consistent  with  previous  GW  studies  using  various  satellite  instruments,  monthly 
mean  AIRS  variance  shows  large  enhancements  over  meridionally  oriented  moun¬ 
tain  ranges  as  well  as  some  islands  at  winter  hemisphere  high  latitudes.  Enhanced 
wave  activities  are  also  found  above  tropical  deep  convective  regions.  GWs  prefer  to 
propagate  westward  above  mountain  ranges,  and  eastward  above  deep  convection. 
15  AIRS  90  field-of-views  (FOVs),  ranging  from  +48°  to  -48°  off  nadir,  can  detect  large- 
amplitude  GWs  with  a  phase  velocity  propagating  preferentially  at  steep  angles  (e.g., 
those  from  orographic  and  convective  sources).  The  annual  cycle  dominates  the  GW 
variances  and  the  preferred  propagation  directions  for  all  latitudes.  Quasi-biennial  os¬ 
cillation  (QBO)  signals  are  also  found  in  the  tropical  lower  stratosphere  despite  their 
20  small  amplitudes. 

From  90  AIRS  FOV  radiance  measurements,  we  are  able  to  clearly  identify  mea¬ 
surement  noises,  high-frequency  internal  GWs,  and  low-frequency  inertia  GWs.  Even 
though  the  vertical  wavelengths  of  inertia  GWs  are  shorter  than  the  thickness  of  in¬ 
strument  weighting  functions,  simulations  support  the  AIRS  sensitivity  to  these  waves. 
25  The  novel  discovery  of  AIRS  capability  of  observing  shallow  inertia  GWs  will  expand 
the  potential  of  satellite  GW  remote  sensing  and  provide  further  constraints  on  the  GW 
drag  parameterization  schemes  in  the  general  circulation  models  (GCMs). 
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1  Introduction 

Gravity  waves  (GWs)  are  known  to  play  a  key  role  in  global  climate  and  weather  dy¬ 
namics  by  transporting  energy  and  momentum  from  the  lower  to  upper  atmosphere, 
transport  which  is  essential  for  determining  the  general  circulation  and  temperature 
5  structure  in  the  stratosphere  and  mesosphere  (e.g.,  mesosphere  wind  reversal,  quasi¬ 
biennial  oscillation,  etc.).  GWs  are  of  particular  importance  for  wave  dynamics  in  the 
summer  hemisphere  in  the  stratosphere  where  planetary  waves  are  weak,  and  are  im¬ 
portant  in  the  entire  mesosphere.  As  an  example,  model  simulation  results  suggest 
that  the  mesosphere  summer  easterlies  would  arrive  one  month  late  without  including 
10  gravity  wave  drag  (GWD)  in  the  model  (Scaife  et  al.,  2002).  GWs  also  affect  weather 
and  chemistry.  For  example,  temperature  fluctuations  associated  with  GWs  can  lead 
to  formation  of  polar  stratospheric  clouds  (PSCs)  in  a  synoptically  warm  condition  and 
susequently  affect  ozone  depletion  (Hamill  and  Toon,  1991). 

The  importance  of  gravity  waves  on  climate  and  weather  became  widely  recognized 
is  in  recent  decades  due  to  progress  in  observational  techniques  and  model  improve¬ 
ments.  Two  issues  remain  as  major  concerns,  however.  Firstly,  most  of  the  GWs 
cannot  be  resolved  in  GCMs  and  hence  need  to  be  parameterized  except  those  with 
long  wavelengths  and  low  frequencies.  These  prescribed  parameters,  formulated  to  be 
dependent  on  detailed  characteristics  of  GWs  and  GW  sources,  however,  are  poorly 
20  constrained  by  observations,  and  are  actually  heavily  “tuned”  in  most  of  the  cases  to 
attain  a  realistic  atmosphere.  For  example,  one  of  the  major  GW  sources,  convection, 
is  often  assumed  to  be  uniformly  distributed  all  over  the  globe  in  GCMs  with  spectral 
GWD  parameterization  schemes,  which  is  obviously  far  from  reality  (Alexander  and 
Rosenlof,  2003).  The  other  important  non-stationary  GW  source  -  jet  imbalances  - 
25  is  not  included  in  most  of  the  GCMs  (Kim  et  al.,  2003).  Secondly,  different  obser¬ 
vational  instruments  can  only  see  partial  gravity  wave  spectra.  Although  solving  this 
jigsaw  puzzle  is  underway  (e.g.,  Alexander  et  al.,  2010),  the  high-frequency  portion  of 
the  spectrum  is  invisible  to  most  of  the  instruments  (Wu  et  al.,  2006;  Alexander  et  al., 
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2010),  and  the  momentum  flux  associated  with  these  waves  remain  uncertain.  Accord¬ 
ing  to  the  gravity  wave  dispersion  relationship,  high-frequency  GWs  tend  to  have  more 
vertical  fluxes  of  the  horizontal  momentum,  and  hence  may  exert  appreciable  impacts 
on  the  mean  flow  with  momentum  deposition.  Therefore,  the  underestimation  of  these 
5  waves  might  result  in  serious  shortfall  for  constraining  the  gravity  wave  drag  (GWD) 
parameterization  schemes  in  current  climate  models. 

Among  various  GW  observations,  satellite  instruments  have  advantages  over  others 
(e.g.,  radiosonde,  radar,  aircraft  measurements)  in  terms  of  high  and  regular  sampling 
density  and  global  coverage.  As  an  example,  the  Aqua  satellite,  which  carries  the 
10  Atmospheric  Infrared  Sounder  (AIRS),  samples  each  2°  x2°  latitude-longitude  gridbox 
at  least  24  times  per  month.  In  addition,  intercomparisons  and  cross-validations  can 
therefore  be  readily  carried  out  among  different  satellite  instruments  on  a  global  ba¬ 
sis.  Based  on  their  viewing  geometry,  satellite  GW  observational  instruments  can  be 
categorized  into  three  groups:  limb,  sub-limb,  and  nadir  sounders.  The  limb  sounders 
is  (e.g.,  LIMS,  CRISTA,  HIRDLS,  GPS)  are  sensitive  to  GWs  with  small  Xz/Xh  ratio,  i.e., 
low-frequency  GWs.  Here  Xz  and  Xx  denote  vertical  and  horizontal  wavelength,  re¬ 
spectively.  The  nadir  sounders  (e.g.,  AIRS,  AMSU-A,  SSMIS)  are  mostly  sensitive 
to  high-frequency  GWs,  and  the  sub-limb  sounders  (e.g.,  MLS)  are  sensitive  to  mid¬ 
frequency  GWs.  More  details  can  be  found  in  Wu  et  al.  (2006). 

20  As  a  nadir  viewing  sounder,  AIRS  has  four  supurb  properties  for  GW  studies.  Firstly, 
its  measured  radiances  are  very  sensitive  to  high-frequency  GWs,  which  are  likely  to 
be  underestimated  so  far  in  both  observations  and  models.  Secondly,  AIRS  radiances 
are  a  direct  measure  of  GW  induced  air  temperature  perturbations,  and  hence  provide 
a  more  accurate  measurement  of  GWs  than  retrieved  temperatures.  Thirdly,  AIRS 
25  has  a  high  horizontal  resolution  (~13km  at  nadir)  which  makes  it  more  attractive  in 
detecting  the  high-frequency,  short  Xh  GWs.  Lastly,  unlike  the  limb  sounders,  AIRS 
sees  GWs  selectively  depending  on  its  viewing  angles.  This  allows  us  to  estimate 
preferred  GW  propagation  directions  from  the  viewing-dependent  variance  difference 
between  the  two  outmost  off-nadir  views. 
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This  paper  is  organized  as  follows.  The  following  section  gives  a  brief  introduction 
of  the  AIRS  instrument,  and  a  description  of  the  methods  we  use  to  retrieve  GW  prop¬ 
erties.  In  Sect.  3,  a  climatology  of  wave  variance  and  the  preferred  zonal  propagation 
direction  will  be  derived  in  terms  of  zonal-mean,  geographical  distribution  and  tempo¬ 
ral  variations.  Interpretation  of  the  results  is  discussed  in  detail  in  Sect.  4,  where  we 
emphasize  on  two  major  gravity  wave  sources  -  topography  and  convection,  and  inter¬ 
annual  variabilities.  Wave  enhancements  in  the  equatorial  lower  stratosphere  and  the 
winter  pole  at  lOhPa  will  be  discussed  in  a  great  length,  and  a  numerical  experiment 
is  performed  to  support  the  speculation  on  the  observed  GWs  being  low-frequency 
waves.  Concluding  remarks  can  be  found  in  Sect.  5. 

For  reader’s  convenience,  we  hereafter  refer  to  “internal  GWs”  as  high-frequency 
(®  »  f),  long  vertical  wavelength  (Az  >  10  km)  and  short  horizontal  wavelength  waves 
{Xh  ~  1 00  km  or  less),  and  use  “inertia  GWs”  to  represent  low-frequency  (®  ~  f),  short 
vertical  wavelength  (Az  <  10  km)  and  long  horizontal  wavelength  {Xh  ~  1000  km)  waves, 
where  ®  and  f  denote  wave  frequency  and  the  Coriolis  parameter.  One  should  keep 
in  mind  that  GWs  in  the  stratified  atmosphere  are  essentially  all  internal  waves. 


2  AIRS  instrument  and  method 

AIRS  is  an  infrared  spectrometer  and  sounder  that  contains  2378  channels  in  3.74- 
4.61  pm,  6.20-8.22  pm  and  8.8-1 5.4  pm  wavebands.  It  makes  a  cross-track  scan  every 
8/3  second  that  includes  90  footprints  on  the  ground,  4  independent  views  over  the 
cold  space,  and  3  views  into  3  different  calibrators  (Aumann  and  Miller,  1 994).  The  scan 
swath  is  ~  1600  km  wide,  reaching  ±48.95°  from  nadir.  The  angle  difference  between 
two  adjacent  footprints,  therefore,  is  about  1.1°,  which  corresponds  to  ~  1 3  km  footprint 
size  at  nadir.  The  satellite  orbits  over  24  h  are  divided  into  240  granules,  each  of  which 
contains  135  scans  (6  min  duration). 

We  use  AIRS  C02  15  pm  radiance  emission  bands  with  wavenumbers  ranging  from 
665  cm-1  to  693 cm-1  in  this  study.  50  channels  are  selected  with  11  distinguished 
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weighting  functions  (WFs)  peaking  between  2  hPa  and  1 00  hPa.  The  channel  numbers 
are  listed  in  Table  A1  in  the  Appendix  A.  Ideally,  only  channels  at  the  crests  or  troughs  of 
the  radiance  spectra,  or  paired  channels  at  two  side-wings  of  peaks  or  troughs  should 
be  used  since  they  are  mostly  stable  with  small  detector  drifts  on  the  spectra,  which 
mainly  occur  due  to  the  expansion/contraction  of  the  glass  filter  depending  on  whether 
it  faces  the  Sun  or  in  darkness  (i.e.,  ascending/descending  of  the  orbits).  However, 
we  are  losing  information  over  most  of  the  middle  to  upper  stratosphere  by  only  se¬ 
lecting  these  channels.  With  this  consideration,  several  channels  on  one  side  wing  of 
peaks/troughs  (bold  numbers  in  Table  A1)  are  also  selected.  Results  from  ascending 
and  descending  orbits  have  been  compared  carefully,  and  no  significant  differences 
have  been  found  whatsoever,  so  this  problem  shouldn’t  cause  major  problems  or  fail¬ 
ures  of  using  the  wing  channels  in  our  research. 

Since  AIRS  scans  are  always  perpendicular  to  the  orbit  track,  the  scan  line  can 
therefore  be  treated  approximately  as  along  the  west-east  direction  for  most  of  the 
time,  except  at  high  latitudes  where  the  scan  is  most  meridional.  We  hence  exclude 
the  data  beyond  ±80°  latitudes. 

The  variance  of  the  radiance  (a2)  has  three  components: 


o2  =  a2  +  o2  .  +  e2 

GW  noise 


(1) 


Normally  other  variance  (e2)  is  much  smaller  than  the  first  terms  on  the  right-hand  of 
Eq.  (1),  and  hence  is  omitted  here.  Evidence  shows  that  AIRS  may  observe  some  tur¬ 
bulence  signals  or  small  vertical  wavelength  inertia  GWs,  which  will  be  discussed  later 
on  in  Sect.  4.3.  However,  these  additional  sources  of  variance  are  highly  uncertain,  and 
we  would  rather  include  them  in  <JqW  than  separating  them  out.  Since  the  instrument 
noise  is  on  the  same  order  with  the  magnitude  of  the  GWs,  it  needs  to  be  carefully 
evaluated  and  subtracted  consequently.  The  instrument  noise  (&„oise)  is  estimated  in 
an  analagous  way  to  Wu  and  Waters  (1996).  Instead  of  using  pre-launch  information, 
we  feel  it  is  more  accurate  to  estimate  the  random  component  of  the  noise  from  the 
real  data  with  the  assumption  that  o20ise  only  depends  on  the  frequency  channel,  and 
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should  not  vary  with  location,  time  and  viewing  angle.  To  extract  the  radiance  noise, 
we  first  apply  a  3-pt  running  smooth  window  to  each  independent  scan  and  compute 
the  variance  as  the  squared  difference  between  the  original  and  smoothed  data  (re¬ 
ferred  as  “3-pt  result”  hereafter).  Then,  the  minimum  variance  from  the  monthly  aver¬ 
aged  2°  x  2°  map  for  each  month  is  determined.  Because  the  3-pt  variance  is  affected 
least  by  GWs,  the  minimum  variance  is  assumed  to  be  very  close  to  the  measure¬ 
ment  noise.  These  minimum  variances  are  found  always  coming  from  the  summer 
hemisphere  high-latitudes.  Finally,  the  mean  is  computed  for  all  the  variances  that  are 
below  the  median  of  the  timeseries,  and  we  attribute  this  mean  as  the  instrument  noise 
CTnoise-  The  ^-pt  variance  is  used  for  the  noise  estimation  since  it  gives  the  smallest 
cut-off  wavelength  and  hence  contains  the  most  high-frequency  information.  Neverthe¬ 
less,  the  3-pt  result  still  contains  a  fair  amount  of  GW  contributions  in  that  the  latitudinal 
and  geographical  distributions  of  the  3-pt  result  look  similar  to  7-pt  result  even  without 
noise  subtraction  (not  shown).  The  noise  estimation  method  used  here  is  confirmed 
to  be  reasonable  and  accurate  by  the  fact  that  variance  reaches  minimum  at  summer 
hemisphere  high  latitudes,  and  is  quite  uniformly  and  stably  distributed  there  through¬ 
out  the  time.  Moreover,  there  are  no  statistically  significant  differences  among  different 
viewing  angles  there. 

We  furthermore  compared  the  noises  estimated  by  the  above  method  with  the  NEdT 
(Noise  Equivalent  Delta  Temperature)  listed  in  AIRS  property  reports  (http://airs.jpl. 
nasa.gov/data_products/algorithms/),  and  as  one  can  easily  tell  from  Table  A1  in  the 
Appendix  A,  they  agree  very  well.  A  linear  increasing  trend  has  been  found  in  the 
estimated  noise  at  various  levels,  which  is  on  the  same  order  of  the  trend  found  in  the 
AIRS  documentation.  Since  this  trend  is  far  from  significant  (regression  coefficients  on 
the  order  of  10“5  K2  month-1),  which  is  probably  due  to  the  run-off  of  the  sensors,  we 
do  not  remove  the  trend  in  this  study. 

As  suggested  by  Wu  and  Waters  (1997)  and  Wu  and  Eckermann  (2008),  the  mini¬ 
mum  detectable  GW  variance  (OqW)  is  o2in  =  \J2{M -2)/A/o20ise,  where  M  is  the  num¬ 
ber  of  points  used  in  this  filter  for  variance  estimation  (M  =  7  in  this  case).  Since  (j20ise 
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is  determined  from  a  monthly  mean  2°  x  2°  map,  a  zonal  average  with  N  points  in  this 
grid  size  can  further  reduce  the  detection  threshold.  The  outmost  M/2  points  are  used 
to  get  an  averaged  variance  value,  which  also  needs  to  be  considered.  The  minimum 
detectable  values  at  various  pressure  levels  are  listed  in  Table  A1  in  the  Appendix  A. 

5  We  cannot  directly  calculate  the  GW  horizontal  propagation  direction  from  radiance 
measurements  of  a  single  scan.  However,  the  preferred  zonal  propagation  direction 
can  be  indirectly  inferred  from  the  difference  of  GW  variances  between  two  different 
FOVs.  This  is  illustrated  in  a  schematic  picture  Fig.  1,  where  a  typical  mountain  GW 
packet  is  generated  downstream  of  the  mountain  under  a  westerly  wind.  Wave  phases, 
io  though  stationary  in  the  Earth  frame,  are  propagating  downward  and  westward  rela¬ 
tive  to  the  westerly  mean  flow,  while  energy  is  propagating  upward  and  westward  in 
this  case.  In  the  west  field-of-view  (FOV)  of  AIRS,  the  wave  amplitudes  are  largely 
smeared  out,  while  the  east  FOV  would  yield  a  much  greater  reading.  Therefore,  nega¬ 
tive  (positive)  value  of  (GWw-GWe)  reflects  a  predominant  westward  (eastward)  zonal 
is  propagation.  Since  multiple  factors  affect  the  magnitude  of  the  difference,  such  as  the 
viewing  angle  relative  to  the  phase  front,  the  original  GW  amplitude,  etc.,  the  absolute 
value  of  the  difference  does  not  have  a  definitive  physical  meaning.  As  long  as  it  is 
statistically  significant,  the  sign  of  the  difference  is  the  variable  that  really  matters.  This 
approach  has  been  formerly  used  on  Aura  MLS  observations  to  study  GW  meridional 
20  propagation  direction,  and  has  turned  to  be  a  very  useful  and  effective  method  (Wu 
and  Eckermann,  2008). 

Before  presenting  the  results,  one  should  note  that  the  majority  of  the  GW  signals  in 
AIRS  comes  from  high-frequency  GWs.  The  visibility  can  be  computed  as  a  function 
of  along-track  wavelength  (Ax)  and  vertical  wavelength  (Az)  by  convolving  the  weighting 
25  function  (WF)  with  a  plane  wave  of  unity  amplitude  (1  K)  (McLandress  et  al.,  2000). 
In  the  along-track  direction,  the  peak  of  the  sensitivity  is  mainly  determined  by  the 
length  of  the  truncation  window,  as  suggested  by  Fig.  2.  The  left  half  of  three  curves, 
which  is  similar  for  the  3-pt,  7-pt  and  15-pt  results,  is  controlled  by  the  width  of  the 
AIRS  FOV.  All  three  rise  sharply  once  Ax  becomes  larger  than  the  footprint  size,  which 
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is  ~  13  km  at  nadir.  The  length  of  the  smoothing  window  decides  the  peak  as  well 
as  the  tail  of  the  filter.  In  the  vertical  direction,  AIRS  can  hardly  detect  GWs  with  Az 
shorter  than  the  full  width  at  half  maximum  (FWHM)  of  the  WFs,  which  is  ~  12  km. 
Flowever,  there  are  some  exceptions,  which  will  be  discussed  later  in  Sect.  4.3.  The 
5  sensitivity  increases  monotonically  with  the  increase  of  the  Az  before  reaching  the  peak 
at  Az  »  Axmax,  where  Axmax  denotes  the  most  sensitive  along-track  wavelength.  This 
can  be  better  understood  from  revisiting  Fig.  1,  where  the  largest  amplitude  would  be 
reached  when  the  GW  fronts  are  parallel  to  the  viewing  angle.  For  a  7-pt  window, 
the  peak  sensitivity  occurs  when  both  Ax  and  Az  equal  to  ~  100  km.  Alexander  and 
io  Barnet  (2007)  studied  multiple  mountain  GW  events  over  the  tip  of  Andes  and  Antarctic 
Peninsula,  and  found  out  the  most  dominant  Ax  is  around  100-1 50  km.  So  the  7- 
pt  window  can  effectively  capture  these  GWs.  With  a  wider  window,  we  are  at  the 
cost  of  possible  inclusion  of  other  waves  (e.g.,  Kelvin  wave),  and  we  may  lose  the 
dominant  GWs  as  suggested  by  Alexander  and  Barnet  (2007).  Nevertheless,  the  major 
15  conclusion’s  robustness  does  not  depend  much  on  window  size.  We  hence  use  7-pt 
results  in  this  paper. 


3  Climatology  of  gravity  wave  variances  and  the  preferred  zonal  propagation 
direction 

The  monthly  mean  latitudinal,  geographical  and  temporal  variations  will  be  presented  in 
20  this  section.  Year  2005  is  used  to  indicate  the  mean  GW  climatology,  which  represents 
a  condition  with  no  strong  sudden  stratospheric  warming  nor  ENSO  events. 

3.1  Latitudinal  distributions 

Figures  3  and  4  plot  the  monthly  averaged  zonal  mean  7-pt  radiance  variance  as  a  func¬ 
tion  of  latitude  and  height  in  January  and  July  of  2005,  respectively  (with  estimated  in- 
25  strumental  noise  deducted),  overplotted  with  UK  Met  Office  monthly  mean  zonal  winds. 
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The  variances  can  be  explained  by  GWs,  since  they  in  general  grow  exponentially  with 
height  at  mid-high  latitudes.  In  the  middle  to  upper  stratosphere  (above  20hPa),  the 
two  largest  variance  centers  colocate  with  the  westerly  and  easterly  jet  cores  in  the 
upper  stratosphere,  and  the  one  in  the  winter  hemisphere  is  apparently  stronger.  This 
5  indicates  that  large  amplitude  GWs  (or  frequently  generated  GWs)  are  from  sources 
beneath  the  two  jet  centers.  The  overall  feature  is  similar  to  previous  GW  studies  using 
Aqua  AMSU-A,  UARS  MLS  and  Aura  MLS  (Wu  et  al.,  2006)  with  comparable  strengths 
in  the  winter  hemisphere,  but  weaker  amplitudes  in  the  summer  hemisphere.  However, 
the  wave  enhancement  at  the  lOhPa  near  the  winter  pole  is  a  new  feature  revealed 
10  by  AIRS  data.  An  interpretation  of  this  feature  is  inertia  GWs,  which  will  be  discussed 
in  Sect.  4.3.  The  GWs  are  almost  tranquil  at  the  summer  polar  region  in  the  AIRS 
observation.  Even  without  the  deduction  of  estimated  noise,  the  variance  there  is  still 
significantly  smaller  than  other  satellite  observations.  AIRS  observations  seem  to  sug¬ 
gest  that  high-frequency  internal  GWs  have  difficulty  propagating  or  being  generated  at 
is  the  mid-high  latitudes  of  the  summer  hemisphere  under  weak  winds.  The  presence  of 
the  zero  wind  line  in  the  summer  hemisphere  is  another  key  reason  whereby  orographic 
GWs  are  mostly  removed. 

In  the  lower  stratosphere,  AIRS  GW  variance  is  generally  larger  in  the  tropics  and 
decreases  with  latitude.  Similar  results  have  been  found  in  some  measurements  that 
20  represent  the  behavior  of  the  inertia  GWs  (e.g.,  Wang  and  Geller,  2003;  Ratnam  et  al., 
2004;  Preusse  et  al.,  2006)  and  middle  frequency  GWs  (e.g.,  Wu  and  Eckermann, 
2008).  AMSU-A,  which  has  the  same  scan  angle  as  AIRS,  and  hence  which  is  sup¬ 
posed  to  be  also  sensitive  to  internal  GWs,  observes  features  opposite  to  AIRS  at  the 
equatorial  lower  stratosphere  where  the  GW  variance  reaches  minimum  (Wu  et  al., 
25  2006).  This  interesting  feature  will  be  discussed  in  details  later  on  in  Sect.  4.3. 

As  stated  in  Sect.  2,  one  can  indirectly  infer  the  GW  preferred  zonal  propagation 
direction  from  the  difference  of  the  GW  variances  between  the  two  outmost  views.  As 
in  the  right  columns  of  Figs.  3  and  4,  the  zonal  mean  zonal  propagation  directions, 
on  a  monthly  scale,  always  tend  to  be  in  the  opposite  direction  to  the  mean  zonal 
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winds,  which  is  in  general  westward  (eastward)  in  the  winter  (summer)  hemisphere. 
This  is  true  especially  for  large  amplitude  GWs  at  high  altitudes,  as  will  be  shown  in 
the  geographical  patterns  in  the  next  section.  Mountain  forcing  is  one  of  the  most 
prominent  GW  sources  in  the  winter  hemisphere.  The  high-frequency  mountain  waves 
generated  in  the  mid-latitude  mountain  areas  can  propagate  exclusively  in  the  relatively 
strong  westerly  wind  to  make  a  zero  ground-based  phase  speed.  The  overall  pattern 
in  Fig.  3c  does  not  strictly  follow  the  mean  zonal  wind  contours,  as  GWs  always  try  to 
drag  the  total  wind,  and  hence  the  meridional  components  of  both  the  wind  and  the  GW 
propagation  need  to  be  taken  into  account.  Note  that  there  is  no  significant  difference 
between  two  outmost  views  at  the  lOhPa  winter  pole  and  tropical  lower  stratosphere, 
which  means  that  the  same  amount  of  eastward  propagating  and  westward  propagat¬ 
ing  GWs  exists  at  those  levels.  Simulation  results  from  Sect.  4.3  confirm  this  inference. 

Very  similar  features  can  also  be  obtained  with  the  3-pt  results,  though  the  ampli¬ 
tudes  are  much  smaller  (not  shown).  Considering  the  fact  that  the  3-pt  running  smooth 
window  is  a  high-pass  filter,  we  can  conclude  that  the  majority  of  GWs  observed  by 
AIRS  are  really  high-frequency,  short  horizontal  wavelength  internal  GWs. 

3.2  Geographical  distributions 

Monthly  mean  geographical  maps  are  plotted  in  Figs.  5  and  6  at  various  levels  to  rep¬ 
resent  typical  situations  of  January  and  July.  Only  the  variance  that  is  larger  than  2crmjn 
is  colored,  where  omin  is  the  minimum  detectable  variance.  Zonal  mean  wind  (black 
solid  contours)  and  wind  speed  at  700  hPa  that  exceeds  10ms_1  (brown  hatched  ar¬ 
eas)  are  overplotted  on  the  variance  maps.  Such  orographic  GWs  have  been  seen 
through  other  satellite  instruments,  such  as  Aura  MLS  (Wu  and  Eckermann,  2008)  and 
AMSU-A  (Wu  et  al.,  2006).  The  GW  amplitudes  in  AIRS,  however,  are  the  strongest 
among  all  even  after  taking  into  account  the  averaging  effect.  This  indicates  that  moun¬ 
tain  GWs  likely  possess  greater  amount  of  energy  toward  the  high-frequency  part  of 
the  frequency  spectrum  (relative  to  the  mean  wind),  and  AIRS  is  particularly  sensitive 
to  them  due  to  its  high  resolution  and  the  steep  viewing  angle.  In  both  hemispheres, 
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orographic  GWs  propagate  poleward  as  they  travel  upward,  in  agreement  with  what 
has  been  found  from  Aura  MLS  observations  (Wu  and  Eckermann,  2008).  The  oro¬ 
graphic  GWs  are  closely  associated  with  the  upper  level  polar  night  jet,  and  hence  the 
location  of  the  enhancements  vary  from  year  to  year  with  the  shift  of  the  polar  night  jet. 
5  Comparing  the  1 0  hPa  map  with  the  2.5  hPa  map,  significant  increase  of  GW  activities 
is  found  in  Alaska  in  January,  2005  (Fig.  5)  and  West  of  Antarctic  Peninsula  in  July, 
2005  (Fig.  6).  They  are  responsible  for  the  peak  at  the  lOhPa  pole  in  the  zonal  mean 
maps  Figs.  3  and  4.  They  are  highly  variable  along  the  longitudes,  and  not  necessarily 
related  with  the  jets.  Interestingly,  the  Brunt-Vaisala  frequency  (A/2)  also  peaks  at  the 
io  same  locations  for  these  two  particular  months  (not  shown).  We  will  discuss  about  this 
feature  later  on  in  Sect.  4.3. 

In  the  subtropics  and  tropics,  large  GW  activities  are  found  in  the  upper  stratosphere 
over  the  deep  convective  regions.  The  deep  convective  regions  are  identified  from 
the  ice  water  content  (IWC)  from  Aura  MLS  (Wu  and  Eckermann,  2008).  In  particu- 
15  lar,  they  are  Western  Pcific  warm  pool  region,  Amazon  rainforest  region,  and  Central 
Africa  rainforest  region  for  NH  winters,  and  Southeastern  US  monsoon  region  and 
India-South  China  monsoon  region  for  SH  winters.  Comparing  with  Aura  MLS,  the 
magnitude  of  these  convectively  generated  GWs  seen  from  AIRS  are  smaller.  This 
is  consistent  with  the  convective  source  spectrum  suggested  by  Beres  et  al.  (2004), 
20  where  source  momentum  fluxes  decrease  rapidly  at  high-frequency  part  of  the  spec¬ 
trum.  The  AIRS  waves  tend  to  propagate  eastward  relative  to  the  mean  wind,  but 
the  differences  between  the  two  outmost  views  have  values  that  are  marginally  above 
the  minimum  detectable  ones,  which  means  both  eastward  and  westward  propagating 
GWs  exist  in  the  upper  stratosphere,  and  the  westward  propagating  ones  are  smaller 
25  than  the  eastward  propagating  ones  with  respect  to  the  integrated  magnitude.  Convec¬ 
tion  is  known  to  generate  a  broad  spectrum  of  GWs  that  propagate  in  both  directions  at 
the  top  of  the  deep  convection  (e.g.,  Alexander  et  al.,  1997;  Beres  et  al.,  2004).  After¬ 
wards,  the  majority  of  the  low-frequency  westward  propagating  GWs  are  wiped  out  by 
the  easterlies  in  the  tropics  and  subtropics,  but  some  of  the  high-frequency  waves  may 
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survive  because  they  generally  have  higher  phase  speeds.  These  convectively  excited 
high-frequency  GWs  do  not  grow  into  a  comparable  strength  with  mountain  GWs  until 
they  reach  the  upper  stratosphere.  In  the  lower  stratosphere,  a  uniform  belt  of  GW 
enhancements  is  found  in  the  tropical  region,  with  some  highlights  in  the  deep  convec¬ 
tive  regions,  which  creates  the  tropical  maximum  in  the  lower  stratosphere  in  the  zonal 
mean  maps.  This  belt  has  been  seen  by  previous  studies  using  various  GW  observa¬ 
tional  instruments  (e.g.,  radiosonde,  GPS,  etc.).  It  apparently  is  not  closely  associated 
with  the  convective  sources.  It  is  very  interesting  that  AIRS  can  see  it  but  AMSU-A,  as 
another  nadir  viewing  sounder,  doesn’t  see  it  at  all.  The  belt  will  be  revisited  later  on  in 
Sect.  4.3,  where  we  claim  the  belt  is  caused  by  the  propagation  of  inertia  GWs  there. 

In  the  high  latitude  summer  hemisphere,  GW  signals  are  barely  detectable.  The  3-pt 
results  without  subtracting  noises  yield  a  quite  uniform  distribution  of  “background-like” 
variance,  which  is  independent  of  location  and  wind.  In  the  NH  summer  (Fig.  6),  it 
is  easy  to  explain  the  NH  serenity  since  the  wind  at  700  hPa  is  too  weak  to  generate 
large  orographic  GWs  (brown  hatched  regions  barely  exist  in  NH  in  Fig.  6).  In  the 
SH  summer,  we  do  see  some  enhancements  over  the  Andes  in  the  lower  stratosphere 
(Fig.  5i),  but  the  waves  hardly  grow  upward  probably  because  of  the  critical  level  filtering 
effect  at  the  zero-wind  line. 

Large  amplitude  mountain  GWs  vary  their  location  year  to  year  (not  shown).  The 
10ms_1  wind  speed  at  700  hPa  roughly  provides  a  good  threshold  for  generating  not- 
icable  topographic  GWs  in  the  stratosphere.  Only  at  places  where  westerly  wind  is 
consistently  strong  from  low  troposphere  to  upper  stratosphere  can  those  GWs  grow 
into  significant  strengths.  Locations  of  high  occurrence  of  convectively  generated  GWs 
are  also  closely  associated  with  the  movement  of  tropical  deep  convective  zones. 

3.3  Temporal  variations 

An  apparent  annual  cycle  with  maximum  in  the  winter  hemisphere  over  mid-high  lat¬ 
itudes  is  found  throughout  the  region  of  interest  in  both  the  variance  and  difference 
timeseries.  We  pick  the  levels  of  2.5,  10,  40  and  80hPa  for  the  timeseries  in  Fig.  7, 
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together  with  the  monthly  zonal  wind  derived  from  UK  Met  Office  dataset.  Since  the 
annual  cycle  is  highly  repeatable,  multi-year  AIRS  variances  are  lumped  together  in 
Fig.  7  as  well  as  the  UKMO  winds.  The  annual  cycle  found  in  AIRS  is  consistent  with 
other  observations  obtained  by  MLS  (Wu  and  Eckermann,  2008),  radiosonde  (Wang 
5  and  Geller,  2003),  GPS  (Tsuda  et  al.,  2000),  etc.  Eckermann  (1995)  attribute  the  an¬ 
nual  cycle  in  both  lidar  and  rocket  sounding  observations  to  seasonal  variations  in  the 
density  stratification  of  the  atmosphere  (i.e. ,  A/2).  Besides  these  major  factors,  source 
properties  can  also  contribute  to  the  variations  in  that  stronger  westerly  winds  dur¬ 
ing  hemispheric  winters  tend  to  create  larger  mountain  GWs,  which  propagate  further 
io  downstream  under  preferable  conditions. 

The  Southern  Hemisphere  differs  from  the  Northern  Hemisphere  in  terms  of  the  lo¬ 
cation  of  the  maximum  variance.  In  the  Northern  Hemisphere,  wave  variance  nearly 
always  peaks  at  the  time  when  jet  is  the  strongest  (December-Feburary),  as  do  the  dif¬ 
ferences.  Interestingly,  in  the  Southern  Hemisphere,  the  variance  peak  occurs  ahead 
is  of  the  peak  of  the  westerly  jet  below  1 0  hPa,  but  the  difference  between  the  two  outmost 
viewpoints  follows  the  variation  of  the  jet  closely,  both  of  which  propagate  northward 
in  time  (May  through  October).  Below  lOhPa,  GW  amplitudes  are  largest  in  July,  but 
the  mean  wind  maximum  comes  in  one  month  later.  A  close  look  at  the  vertical  config¬ 
uration  of  the  zonal  wind  reveals  that  tropospheric  westerlies  in  the  SH  high  latitudes 
20  actually  peak  in  July  (not  shown).  Another  plausible  explanation  is  based  on  the  fact 
that  the  coastline  of  Antarctica  is  approximately  parallel  to  the  latitude  lines.  Therefore, 
meridional  winds  are  more  important  than  zonal  winds  there  in  generating  large  ampli¬ 
tude  waves.  Actually,  the  meridional  wind  contours  follow  the  change  of  GW  variance 
more  closely  than  zonal  wind  contours  (not  shown).  The  Andes  seem  to  be  vital  at 
25  high  levels  (above  lOhPa)  at  the  intensification  stage  of  the  jet  speed,  while  Antarctic 
Peninsula  is  more  important  in  lower  levels  and  seems  to  be  a  major  factor  causing 
the  northward  movement  of  the  center  of  polar  night  jet.  This  phenomenon  is  also  ev¬ 
ident  in  Fig.  6.  Note  that  the  polar  night  jet  in  SH  winter  (e.g.,  Fig.  3)  curves  poleward 
(equatorward)  below  (above)  lOhPa.  Antarctic  Peninsula  and  the  Andes  perhaps  play 
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different  roles  at  different  levels  in  the  formation  of  this  shape.  In  the  NH  winter,  there 
is  no  such  northward  curvature  in  either  the  jet  center  or  the  center  of  wave  variance 
maximums.  This  result  again  suggests  the  importance  of  mountain  GWD  in  shaping 
and  regulating  the  polar  night  jets. 

5  Strong  interannual  variability  exists  in  NH  high  latitudes  (not  shown).  Polar  night  jet 
splits  during  the  winters  of  2004,  2006  and  2009  at  the  polar  region,  and  GW  variance 
reduces  a  significant  amount  simultaneously.  These  are  the  years  of  strong  sudden 
stratospheric  warming  (SSW)  events  (Lee  et  al.,  2009).  Lee  et  al.  (2011)  observed  en¬ 
hanced  planetary  wave  activities  during  the  SSW  events  from  Aura  MLS  observations, 
io  and  they  believe  the  planetary  waves  serve  as  an  efficient  filter  to  prevent  upward  prop¬ 
agation  of  GWs.  More  explanations  can  be  found  in  Dunkerton  and  Butchart  (1984). 

The  annual  cycle  also  dominates  the  subtropical  stratosphere,  but  completely  oppo¬ 
site  phases  are  found  between  the  upper  stratosphere  and  lower  stratosphere.  In  the 
upper  stratosphere  (first  row  of  Fig.  7),  GW  variance  as  well  as  difference  peak  during 
is  hemispheric  summers.  Zonal  mean  maps  (Figs.  5  and  6)  show  that  the  convectively 
generated  GWs  grow  significantly  in  the  upper  stratosphere,  which  collocate  with  the 
easterly  jet  centers  as  well.  Therefore,  this  summer  peak  could  be  largely  explained 
by  the  convective  source  properties.  In  the  lower  stratosphere,  winter  maximum  dom¬ 
inates,  and  this  should  be  again  attributed  to  the  variations  of  density  stratification 
20  similar  to  that  of  the  mid-high  latitudes  (Eckermann,  1 995). 

The  tropical  lower  stratosphere  has  an  annual  cycle  with  NH  winter  maximum,  and 
this  largely  agrees  with  the  fact  that  deep  convections  in  the  SH  summers  are  in  general 
stronger  (Liu  et  al.,  2007).  QBO  signals  are  relatively  weak  but  still  significant  compared 
with  the  annual  cycle,  showing  very  interesting  features  that  will  be  discussed  in  the 
25  next  section.  In  the  middle  to  upper  stratosphere,  a  prevailing  semi-annual  cycle  is 
found  (not  shown),  which  we  believe  is  a  combined  effect  of  background  winds  (semi¬ 
annual  oscillations  there)  and  the  source  (convection  crossing  the  equator  twice  per 
year  at  equinoxes). 
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4  Discussions 

Some  new  features  from  AIRS  GW  observations  are  worth  further  discussions.  In  this 
section,  we  will  provide  in-depth  analyses  on  GW  sources,  interannual  variations,  and 
enhancements  from  inertia  GWs. 

5  4.1  Orographic  and  convective  GW  sources 

Topographic  GWs  have  been  extensively  studied  during  the  last  decades,  and  relatively 
mature  parameterization  schemes  of  the  mountain  GWDs  have  already  been  employed 
widely  in  state-of-art  GCMs.  These  schemes  usually  take  into  account  factors  such  as 
mountain  height,  surface  wind  speed,  atmospheric  stability  and  the  alignment  of  the 
io  mountain  ridges  (McFarlane,  1987;  Bachmeister,  1993).  The  non-stationary  convec¬ 
tive  GWs  received  less  attention  compared  to  the  orographic  sources  until  recent  years 
(Song  and  Chun,  2005;  Beres  et  al.,  2005).  As  to  the  parameterization  of  these  GWs, 
some  treat  them  as  “moving  mountains”,  and  some  think  they  create  much  broader 
spectra  (e.g.,  Song  and  Chun,  2005).  Both  of  the  orographic  and  convective  GWD  pa- 
15  rameterizations  have  free  parameters  to  represent  the  wave  intermittency,  which  is,  by 
far,  the  poorest  constrained  parameter.  Other  “tunable”  parameters  include  horizontal 
wavelength,  source  momentum  spectrum,  etc.  AIRS,  as  one  powerful  GW  observa¬ 
tion  instrument,  can  certainly  improve  our  global  understanding  on  the  “observational 
constraints”  of  the  parameterization  schemes.  For  example,  from  the  geographical  dis- 
20  tribution  that  has  been  discussed  in  Sect.  3.2,  we  can  tell  that  mountain  GWs  tend 
to  have  larger  energy  variance  toward  the  high-frequency  part  of  the  energy  spectrum, 
while  it  seems  to  be  the  other  way  around  for  the  convective  GWs.  Besides,  some  inter¬ 
esting  features  are  found  in  AIRS,  which  are  quite  different  from  what’s  been  observed 
before  by  other  instruments. 

25  To  look  into  more  details,  we  pick  two  orographic  GW  cases  and  one  convective 
GW  case  that  are  evident  in  January  2005,  which  occur  at  the  Rockies,  the  warmpool 
region,  and  Iceland.  The  monthly  mean  FOV-dependent  GW  variance  within  the  white 
grid  box  is  plotted  at  2.5  hPa  and  80  hPa  for  both  January,  2005  and  2008  in  Fig.  8.  At 
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first  glance,  the  linear  regression  line  confirms  our  idea  that  orographic  GWs  in  general 
propagate  westward  while  tropical  convective  GWs  propagate  eastward,  both  relatively 
to  the  local  mean  zonal  wind.  Ideally,  if  a  monochromatic  wave  exists,  we  should 
see  a  smooth  curve  with  a  single  peak  or  no  peak  depending  on  whether  the  upward 
5  propagation  angle  is  smaller  or  larger  than  48°.  However,  multiple  peaks  stand  out  in 
almost  all  curves  in  Fig.8d-p,  indicating  a  spectrum  of  GWs  are  generated  down  below. 
In  the  lower  stratosphere  (Fig.  8k-p),  convectively  excited  GWs  may  have  sporadic 
spikes  (Fig.  8i,o)  while  orographic  GWs  can  generate  a  rather  broad  spectrum  (Fig.  8n). 
The  peaks  of  the  variance  in  certain  cases  are  close  to  the  nadir  (e.g.,  Fig.  8o,p), 
10  which  means  those  Doppler-shifted  GWs  propagate  almost  vertically.  By  comparing 
the  2.5  hPa  FOV  curves  with  the  ones  at  80hPa,  we  may  conclude  that  background 
wind  plays  a  leading  role  in  shaping  the  GW  spectrum  during  the  upward  propagation. 
However,  three  preferred  angles  at  FOV  No.  25,  No.  0  and  No.  20  persist  from  lower  to 
upper  altitudes  over  the  Rockies  during  January  2005,  and  the  FOV  No.  0  and  No.  20 
15  also  have  peaks  of  GW  variance  at  2.5  hPa  during  January  2008,  even  though  no  big 
wave  event  occurred  over  the  Rockies  in  January  2008.  Similar  phenomena  happen 
over  Iceland  (Fig.  8f,j),  and  the  Andes  (not  shown).  These  preferred  directions  should 
be  associated  with  the  detailed  structures  of  the  topography,  and  background  wind  is 
probably  another  candidate  that  prefers  to  select  certain  angles. 

20  GWs  seen  by  AIRS  provide  at  least  two  new  discoveries  about  the  orographic  and 
convective  GWs.  For  orographic  GWs,  vertical  wavelength  is  already  large  at  the  wave 
generation  level  compared  with  convective  GWs,  and  certain  angles  are  favored  for 
generating  large  amplitude  GWs,  which  vary  among  different  mountain  ranges  and 
vertical  wind  structures  as  well.  The  ratio  between  peaks  and  troughs  as  well  as  max- 
25  imum  amplitude  occurring  FOV  angles  can  be  used  to  infer  wave  intermittency,  which 
is  an  important  wave  parameter  in  GCMs  that  lacks  observational  constraints.  It  is  not 
clear  why  certain  angles  are  preferred.  The  reason  why  large  amplitude  GWs  can  be 
seen  in  the  nadir  view  is  another  interesting  point  and  beyond  the  scope  of  this  study. 
They  remain  as  potential  topics  for  future  investigations. 
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Convectively  generated  GWs,  as  suggested  here  in  AIRS  observations,  have  much 
smaller  wave  scales  in  terms  of  both  vertical  and  horizontal  wavelengths.  Therefore, 
it  is  expected  to  have  fair  amount  of  momentum  flux  in  the  high-frequency  end  of  the 
momentum  spectrum,  which  also  varies  quite  a  lot  based  on  sporadic  spikes  present 
5  in  AIRS.  The  current  assumptions  of  source  spectrum  (e.g.,  Song  and  Chun,  2005; 
Beres  et  al.,  2004)  for  convective  GWs  is,  at  least,  not  accurate  enough  to  represent  all 
situations.  Note  that  these  curves  are  all  derived  from  monthly  means.  Spikes  in  both 
mountain  and  convection  cases  may  simply  come  from  different  events.  Single  case 
studies  need  to  be  carried  out  in  the  future  with  sophisticated  model  runs  to  give  more 
io  definitive  conclusions. 

4.2  Interannual  variations  at  the  equatorial  region 

The  quasi-biennual  oscillation  (QBO)  acts  as  an  efficient  GW  filter  in  the  tropical  lower- 
middle  stratosphere,  and  the  interations  between  GWs  and  background  winds  are,  in 
turn,  an  essential  ingredient  in  the  QBO.  Different  GWs  contribute  differently  to  driving 
is  the  QBO.  Kawatani  et  al.  (2010a,b)  were  able  to  simulate  a  realistic  QBO  with  a  high- 
resolution  climate  model  without  parameterizing  GWD.  Their  GCM  can  resolve  part  of 
the  inertia  GW  spectrum  (limited  by  the  model’s  vertical  grid  resolution),  and  the  sim¬ 
ulated  QBO  has  a  shorter  period  than  that  of  the  real  atmosphere,  but  with  realistic 
strength  and  structure.  Their  work  suggests  the  importance  of  inertia  GWs  in  gen- 
20  erating  the  QBO.  Some  other  model  studies  suggest  that  the  momentum  deposition 
associated  with  the  filtered  GWs  accounts  for  at  least  half  of  the  total  momentum  flux 
in  driving  the  QBO  (Giorgetta  et  al.,  2002;  Kawatani  et  al.,  2010a).  QBO  signals  have 
also  been  seen  in  various  GW  observations,  such  as  radiosonde  (Vincent  and  Alexan¬ 
der,  2000;  Wang  and  Geller,  2003),  GPS  (Torre  et  al.,  2006)  and  Aura  MLS  (Wu  and 
25  Eckermann,  2008).  These  instruments  are  ideal  for  seeing  low  to  middle  frequency 
GWs. 
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AIRS  GWs  contain  a  strong  annual  cycle  in  the  tropical  lower  stratosphere,  as  shown 
in  Fig.  7.  On  top  of  the  annual  cycle,  a  clear  QBO  signal  stands  out.  After  removing  the 
annual  cycle,  the  linear  trend,  and  sub-seasonal  oscillations,  we  present  the  timeseries 
of  the  GW  variance/difference  anomalies  at  the  equator  in  Fig.  9.  Apparent  GW  vari- 
5  ance  enhancements  occur  along  with  the  descent  of  the  QBO  easterly  phase  (defined 
as  when  the  monthly  zonal  wind  is  easterly),  and  a  reduction  of  GW  activity  happens 
right  before  and  lasts  through  the  whole  period  of  the  descending  of  the  westerly  phase 
(Fig.  9a).  Year  2003-2004  is  somewhat  different  from  other  years,  when  the  enhance¬ 
ment  lasts  for  a  very  short  time,  and  a  strong  negative  anomaly  continues  throughout 
io  the  descending  of  the  QBO  easterly  shear.  The  GW  differences  also  show  preferred 
zonal  propagation  directions  during  different  QBO  phases.  Although  it  is  more  am¬ 
biguous  in  Fig.  9b,  in  general  the  GWs  that  can  survive  the  wind  perfer  to  propagate 
eastward  (westward)  in  the  presence  of  QBO  easterly  (westerly)  shears.  Again,  year 
2003-2004  is  an  exception. 

is  The  phenomena  of  enhanced  (reduced)  GW  variance  in  the  QBO  easterly  (westerly) 
phase  can  be  explained  as  follows.  Deep  organized  convection  is  the  most  dominant 
source  at  the  tropics  that  is  responsible  for  GWs  we  see  in  AIRS,  as  discussed  in 
the  last  section.  Since  most  of  the  deep  convective  events  in  the  equatorial  region 
move  eastward  (Wang  and  Rui,  1990;  e.g.,  convections  embedded  in  Madden-Julian 
20  Oscillations  during  boreal  winters  and  Intra-Seasonal  Oscillations  during  boreal  sum¬ 
mers),  they  tend  to  generate  eastward  propagating  GWs.  In  the  QBO  easterly  phase, 
the  eastward  propagating  waves  can  survive,  while  they  are  mostly  filtered  out  in  the 
QBO  westerly  phase.  This  leads  to  the  observed  appearance  that  more  GWs  sur¬ 
vive  in  QBO  easterly  phase.  This  result  is  consistent  with  what  has  been  observed  by 
25  Aura  MLS  (Wu  and  Eckermann,  2006),  where  suppressions  (enhancements)  of  wave 
activities  are  observed  in  the  QBO  westerly  (easterly)  phase  in  the  equatorial  lower 
stratosphere.  It  suggests  that  the  GWs  observed  by  AIRS  may  play  a  more  important 
role  for  the  descent  of  QBO  westerly  phase  than  that  of  the  easterly  phase  as  more 
AIRS  GWs  are  removed  and  hence  deposit  their  momentum  fluxes  in  the  QBO  westerly 
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phase.  However,  the  AIRS  variance  explained  by  the  QBO  is  about  10  times  smaller 
than  that  explained  by  the  annual  cycle,  which  means  the  GW  variance  is  only  slightly 
modified  by  the  QBO  rather  than  playing  a  dominant  role  on  the  formation/propagation 
of  QBO  phases.  This  is  expected  since  GWs  seen  by  AIRS  are  mostly  high  frequency 
5  waves  that  are  usually  with  fast  vertical  group  velocity.  They  usually  break  in  the  upper 
stratosphere  and  above,  and  significantly  modulate  the  mean  wind  there  (Alexander 
and  Fritts,  2003). 

Some  of  the  previous  studies  using  radiosonde  and  GPS  also  found  suppression  of 
GWs  in  the  QBO  westerly  phase  (e.g.,  Sato  and  Dunkerton,  1997;  Wang  and  Geller, 
io  2003;  Torre  et  al.,  2006),  but  they  suggest  that  GWs  are  significantly  enhanced  in  the 
QBO’s  descending  westerly  shear  zone.  Since  Kelvin  waves  are  unavoidably  blended 
in  the  measurements  from  radiosonde  and  GPS,  they  are  a  major  contributor,  and  help 
to  bring  the  westerly  phase  down  because  their  momentum  flux  is  eastward.  MLS 
scans  along-track  (almost  meridionally  over  most  of  the  latitudes),  and  Kelvin  waves 
is  are  truncated  by  a  technique  with  short  horizontal  wavelength  cut-off.  7-pt  running 
smooth  window  we  applied  on  AIRS  measurements  produces  a  cut-off  wavelength  at 
~  105  km,  which  can  also  effectively  filter  out  Kelvin  waves.  Another  possible  reason 
lies  in  the  fact  that  radiosonde  and  GPS  are  likely  observing  different  parts  of  the  GW 
spectrum  from  AIRS.  Vincent  and  Alexander  (2000)  found  increased  GW  activities  dur- 
20  ing  QBO  easterly  phase  through  radiosonde  data  over  a  tropical  island,  which  is  a  bit 
different  from  other  radiosonde  observations  but  is  consistent  with  what  we  observe  in 
AIRS.  The  reason  for  the  strong  negative  anomaly  during  year  2003-2004  remains  un¬ 
clear.  It  might  be  associated  with  other  interannual  variabilities  such  as  El  Nino  Sothern 
Oscillation  (ENSO). 

25  Theoretical  work  predicts  both  westward  and  eastward  propagating  GWs  should  both 
contribute  to  the  QBO  formation  (Baldwin  et  al.,  2001).  This  is  further  confirmed  by 
simulation  results  from  high-resolution  GCMs  (e.g.,  Sato  et  al.,  1999;  Watakani  et  al., 
2010a,b).  However,  the  GWs  with  eastward  phase  speed  are  hardly  observed  through 
various  GW  measurement  techniques  (e.g.,  Sato  et  al.,  1997;  Wheeler  et  al.,  2000)  or 
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it  is  difficult  to  separate  these  GWs  from  massive  Kelvin  wave  signals  as  stated  in  the 
last  paragraph.  AIRS  observes  the  tilt  of  wave  phases  in  the  zonal-height  plane,  from 
which  the  zonal  direction  of  the  intrinsic  phase  propagation  can  be  inferred,  assuming 
that  the  group  propagation  is  upward.  The  inferred  preferred  zonal  wave  propagation 
direction  from  AIRS  at  the  equator  show  both  eastward  and  westward  propagation  sig¬ 
nals  with  about  equal  probability  to  occur  (Fig.  9b).  In  particular,  since  GWs  with  phase 
speed  of  the  same  sign  with  the  zonal  wind  tend  to  meet  the  “critical  level”  more  easily, 
westward  (eastward)  propagating  GWs  are  more  abundant  in  QBO  westerly  (easterly) 
shear.  This  is  the  first  time  the  eastward  propagating  GWs  are  seen  in  satellite  GW 
observations.  It  indicates  the  powerfulness  and  uniqueness  of  satellite  instruments  in 
obtaining  the  GW  information. 

4.3  Wave  enhancement  in  the  equatorial  lower  stratosphere  and  near  the  winter 
pole  at  lOhPa 

Mountain  GWs  are  excessively  strong  in  the  AIRS  observations,  which  overwhelm 
other  GW  components.  In  the  lower  stratosphere,  however,  GWs  are  stronger  in  the 
equatorial  region.  These  GWs  are  not  closely  associated  with  deep  convetion,  as  the 
variance  forms  a  uniform  belt  (Fig.  5h,i,  Fig.  6h,i,  and  bottom  row  of  Fig.  7).  By  revisit¬ 
ing  Fig.  3,  we  can  find  that  the  amplitude  of  these  GWs  maximizes  at  the  tropopause 
(~  lOOhPa),  decreases  dramatically  right  above,  and  increases  again  above  80hPa. 
These  features  are  in  good  agreement  with  previous  studies  using  radiosonde  (e.g., 
Wang  and  Geller,  2003),  GPS  (e.g.,  Ratnam  et  al.,  2004),  Aura  MLS  (Wu  and  Ecker- 
mann,  2008),  CRISTA  and  SABER  (Preusse  et  al.,  2006),  and  high-resolution  GCM 
simulations  (e.g.,  Sato  et  al.,  1999),  but  not  shown  in  UARS  MLS  (Wu  and  Waters, 
1996)  nor  Aura  AMSU-A  (Wu,  2004).  Meanwhile,  there  is  also  a  peak  at  the  winter 
pole  at  1 0  hPa,  as  shown  in  Figs.  3  and  4.  This  is  the  first  observation  of  such  a  1 0  hPa 
winter-pole-enhancement.  Since  this  feature  is  not  uniformly  distributed  among  the 
high-latitude  hemisphere  winters  (Figs.  5  and  6),  and  the  locations  vary  from  year  to 
year  (not  shown),  it  is  believed  not  purely  caused  by  underestimation  of  the  instrument 
noise,  but  rather  a  real  atmospheric  phenomenon. 
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A  careful  evaluation  of  FOV-dependent  variance  profile  (“FOV  curve”  hereafter)  sug¬ 
gests  that  the  lOhPa  variance  is  something  other  than  high-frequency  internal  GWs 
and  instrumental  noise.  The  zonal  mean  FOV  curves  derived  at  60°  N,  0°  and  60°  S 
for  January,  2005  at  80hPa  are  plotted  in  Fig.  10.  At  60°  N,  as  we  would  expect  to 
5  see  from  topographically  generated  internal  GWs  in  Northern  Plemisphere  winter,  the 
variance  is  the  largest  at  the  eastmost  view.  At  60°  S,  the  FOV  curves  are  basically 
flat  and  small,  close  to  the  estimated  noise  levei.  At  the  equator,  the  FOV  curve  bends 
downward,  with  maximum  value  occurring  at  the  nadir  view.  This  is  also  true  at  the  two 
flanks  of  the  polar  night  jet  at  middle  to  upper  stratosphere  (not  shown). 
io  Turbulence  can  cause  “bell”  shape  FOV  curves  seen  in  Fig.  10b  due  to  the  fact 
that  the  power  index  for  2-D  turbulence  is  more  negative  at  the  off-nadir  view  than 
the  nadir  view  and  hence  the  area  integral  is  smaller  for  the  off-nadir  view  (Gruninger 
et  al.,1998).  Since  shear  instabilities  easily  occur  at  the  two  flanks  of  the  polar  night  jet, 
turbulence  can  be  particularly  strong  there.  Plowever,  we  cannot  explain  why  turbulence 
is  is  particularly  strong  in  the  equatorial  lower  stratosphere  and  the  winter  pole  at  1 0  hPa 
with  this  theory.  Besides,  according  to  Aumann  and  Miller  (1994),  the  most  off-nadir 
view  has  a  FOV  swath  area  4  times  larger  than  that  of  the  nadir  view,  which  could 
result  in  a  deeply  downward  curved  FOV-dependent  variance  profile.  Even  with  noise 
being  already  subtracted,  the  FOV  curve  in  Fig.  10b  has  a  difference  between  off-nadir 
20  and  nadir  views  of  only  7%,  which  is  way  too  small  to  be  purely  interpreted  as  the 
turbulence. 

Another  explanation  of  the  downward  curvature  of  FOV  curves  in  Fig.  10b  invokes 
low-frequency  inertia  GWs,  and  we  will  show  that  AIRS  can  in  fact  observe  them. 
Alexander  et  al.  (2002)  proposed  the  idea  that  slow  GWs  have  larger  probabilities  to 
25  be  observed  than  fast  waves.  Since  the  vertical  group  velocity  is  approximately  propor¬ 
tional  to  square  of  the  wave  intrinsic  frequency,  AIRS  has  a  larger  chance  to  observe 
inertia  GWs  in  the  lower  stratosphere,  and  the  overall  enhancement  of  wave  activity 
observed  in  the  tropics  is  largely  due  to  the  reason  that  these  inertia  GWs  propagate 
slantwise  and  spread  out  from  the  source  consequently.  Since  the  GW  enhancement 


ACPD 

11, 11691-11738,  2011 


Gravity  wave 
information  derived 
from  AIRS  radiances 

J.  Gong  et  al. 


Title  Page 


Abstract 


Introduction 


Conclusions 


References 


Tables 


Figures 


Back 


Close 


Full  Screen  /  Esc 


Printer-friendly  Version 


Interactive  Discussion 


11712 


in  the  equatorial  lower  stratosphere  observed  by  AIRS  is  similar  to  GWs  revealed  by  ra¬ 
diosonde  and  GPS  as  well  as  simulation  results  in  high-resolution  GCMs  (e.g.,  Fig.  8c 
in  Sato  et  al.,  1 999),  we  have  more  confidence  to  claim  that  AIRS  observes  inertia  GW 
signals. 

5  Here  we  conduct  a  sensitivity  simulation  to  show  how  AIRS  variance  has  a  good  sen¬ 
sitivity  to  inertia  GWs  in  the  lower  stratosphere.  AIRS’s  vertical  WFs  have  a  thickness 
of  about  12  km,  which  leads  to  a  strong  response  to  long  vertical  wavelength  waves 
(Alexander  and  Barnet,  2007).  The  WF  is  practically  a  delta  function  with  respect  to 
these  long  vertical  wavelength  GWs  and  therefore  yield  GW  variance  at  the  level  where 
io  the  WF  peaks.  For  waves  with  vertical  wavelength  less  than  12  km,  the  convolution  of 
AIRS  WF  with  waves  of  a  constant  amplitude  will  yield  a  very  small  variance.  However, 
if  the  wave  amplitude  varies  with  height,  the  GWs  at  the  edge  of  the  WF  cannot  be 
completely  cancelled,  and  will  yield  some  response  in  AIRS  variance.  Figure  1 1  a  gives 
an  ideal  WF  at  nadir  view  at  80hPa  to  represent  mean  WFs  of  multiple  AIRS  chan- 
15  nels  we  select  which  all  peak  at  80  hPa  (Appendix  A).  The  WF  at  the  outmost  off-nadir 
view  covers  an  area  4  times  broader  than  that  of  the  nadir  WF  (see  Appendix  B  for  de¬ 
tails).  Figure  1 1  b  shows  the  GW  that  is  a  result  of  the  original  imported  GW  convolved 
with  the  WF.  The  original  GW  has  an  amplitude  varying  with  A/2,  and  the  maximum 
amplitude  Amax  is  5  K  (see  Appendix  B  for  details). 

20  The  final  goal  of  this  experiment  is  to  try  to  find  the  best  wave  parameters  that  give 
the  convolved  variance  at  the  nadir  and  the  variance  difference  between  the  nadir  and 
outmost  views  consistent  with  the  observations  (black  line  in  Fig.  10b).  In  order  to 
compute  the  variance,  we  make  a  horizontal  phase  shift  of  mz  forward  of  the  convolved 
wave,  as  shown  in  shades  of  Fig.  1 1  b.  This  procedure  is  to  account  for  the  7-pt  smooth- 
25  ing  window  we  applied  to  AIRS  data.  The  difference  between  the  two  convolved  waves 
are  then  calculated  and  compared  with  the  observation  (Fig.  10b)  at  the  nadir  view. 
With  a  fixed  Xx  of  700  km,  7-pt  smoothing  window  length  (~  105  km)  corresponds  to 
0.1 57T  phase  shift.  We  can  therefore  plot  out  the  expected  AIRS  variance  as  a  func¬ 
tion  of  vertical  wavelength  Xz  in  Fig.  12a.  The  same  procedure  can  be  applied  to  the 
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extreme  off-nadir  point  to  obtain  a  reading  of  wave  variance  there,  and  the  difference 
between  the  nadir  and  outmost  views  can  also  be  computed  as  shown  in  Fig.  12b.  The 
two  curves  meet  with  the  observed  values  both  at  A,  =  6  km.  We  tried  other  Av  values 
and  it  turns  out  Ax  =  700  km,  Az  =  5  km  and  Amax  =  5  K  gives  the  best  agreements  simul- 
5  taneously.  Typical  horizontal  and  vertical  wavelengths  for  equatorial  inertia  GWs  are 
around  1000  km  and  5  km,  respectively  (Wang  and  Geller,  2003),  and  the  peak  ampli¬ 
tude  is  around  5K  (Sato  et  al.,  1999),  which  is  in  good  agreement  with  our  simulation 
result.  The  inertia  GW  variance  can  also  explain  no  obvious  propagation  direction  in 
the  difference  between  eastmost  and  westmost  views,  as  the  radiance  response  in- 
io  volves  no  propagation  direction  information  whatsoever  in  this  case.  The  layer  near 
80hPa  (Fig.  11a)  is  less  affected  by  the  QBO  wind  than  the  layers  above  and  hence 
we  still  do  not  see  a  dominant  QBO  signal  even  these  waves  are  believed  to  be  dom¬ 
inantly  inertia  GWs.  The  simulation  with  a  realistic  AIRS  WF  supports  the  idea  that 
the  enhanced  GW  variance  at  lOhPa  winter  pole  is  also  low-frequency  inertia-gravity 
is  waves.  See  Appendix  B  for  details  about  this  experiment. 

The  inertia  GWs  at  the  equatorial  lower  stratosphere  cannot  propagate  off-equator 
further  to  a  high  altitude  since  their  paths  are  too  slantwise,  and  they  easily  meet  the 
lower  limit  of  frequency,  which  is  the  inertia  frequency  f,  as  f  increases  with  latitude 
(Sato  et  al.,  1999).  Hence  the  tropical  belt  of  GW  enhancement  disappears  at  high  alti- 
20  tudes  (20  hPa  and  above),  and  only  those  high-frequency  internal  GWs  that  are  closely 
related  with  convective  sources  can  survive.  Choi  et  al.  (2011)  applied  3-D  AIRS  WFs 
to  their  parameterized  convective  GWs,  and  found  out  that  merely  any  GWs  can  be 
seen  in  the  equatorial  lower  stratosphere  since  the  parameterized  GWs  are  too  small 
in  terms  of  the  vertical  wavelength  and  horizontal  wavelength  as  well.  Propagating 
25  upward,  these  small  GWs  are  Doppler-shifted  toward  longer  vertical  wavelength,  and 
become  detectable  by  AIRS,  which  agree  very  well  with  what  we  see  at  middle  to  upper 
stratosphere. 

A  peak  at  1 0  hPa  at  the  winter  pole  appears  in  the  A/2  field  in  Fig.  1  of  the  Appendix  B. 
Therefore,  as  in  the  tropical  lower  stratosphere,  inertia  GWs  should  grow  again  at  those 
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levels,  but  with  a  short  penetration.  Convolving  with  the  weighting  function,  and  the 
GWs  result  in  another  enhancement  in  AIRS  variance.  Due  to  the  sensitivity  limitations, 
radiosonde  and  GPS  sensors  cannot  observe  those  GWs  at  this  level.  UARS  MLS  and 
Aura  AMSU-A  cannot  see  these  features  probably  because  of  the  coarse  horizontal 
resolutions  (Wu  and  Eckermann,  2008).  Aura  MLS  also  observes  similar  equatorial 
enhancement  at  lower  levels,  but  not  at  the  lOhPa  winter  pole,  which  remains  to  be 
investigated  in  the  future  (Wu  and  Eckermann,  2008). 

After  all,  the  AIRS  GW  variance  is  believed  to  be  composed  of  large  amplitude  in¬ 
ternal  GWs  that  are  closely  related  with  the  local  sources  (e.g.,  orography,  deep  con¬ 
vections),  and  small  amplitude  “background-like”  inertia  GWs  that  are  determined  by 
the  shape  of  the  WFs  and  the  shape  of  the  original  amplitude  spectra.  It  is  shown 
that  AIRS  can  observe  the  inertia  GWs  even  though  their  vertical  wavelengths  are  less 
than  AIRS  WF  thickness.  This  interpretation  might  be  applicable  to  explaining  signals 
discovered  in  other  instrument  measurements.  The  two  groups  of  GWs  are  potentially 
separable  by  regressing  the  variance  to  A/2.  Turbulence  is  another  possibility,  but  the 
strength  and  mechanism  remain  as  a  question  for  future  investigation. 

c/5' 

o 

c 
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5  Conclusions 

In  this  study,  GW  induced  temperature  variance  is  succesfully  extracted  from  the  bright¬ 
ness  temperature  variance  observed  by  Aqua  AIRS  by  carefully  removing  instrumental 
noises  at  various  pressure  levels.  A  total  of  50  AIRS  channels  is  used  with  WFs  peak¬ 
ing  at  pressure  levels  between  2hPa  and  lOOhPa,  covering  GWs  in  the  entire  strato¬ 
sphere.  Because  of  the  nadir  viewing  geometry  and  small  footprint  size,  AIRS  radiance 
variance  is  mostly  sensitive  to  high-frequency  internal  GWs  with  comparable  vertical 
and  horizontal  wavelengths  with  wave  fronts  in  parallel  to  the  scanning  angle.  Large 
smearing  occurs  when  the  wave  front  is  across  the  AIRS  scan,  leading  to  possibility 
of  inferring  the  preferred  GW  propagation  direction  using  variance  difference  between 
two  outmost  views.  Considering  the  fact  that  the  scan  is  perpendicular  to  the  orbit 
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paths  and  AIRS  is  in  a  polar  orbit,  the  inferred  wave  propagation  direction  is  mostly  in 
the  zonal  direction.  Climatology  of  the  GW  variance  and  zonal  preferred  propagation 
direction  are  documented  in  this  study. 

Comparisons  have  been  carried  out  throughout  the  paper  between  AIRS  GW  vari- 
5  ance  and  other  measurements  including  satellite  and  non-satellite  instruments  (e.g., 
Aura  MLS,  Aqua  AMSU-A,  GPS,  radiosonde,  GCM  simulations,  etc.).  The  GW  vari¬ 
ance  seen  by  AIRS  in  general  agrees  with  previous  studies  in  terms  of  the  distributions 
and  variations  of  maxima.  The  AIRS  GW  peaks  are  generally  found  to  be  above  high 
meridionally-oriented  mountain  ranges,  as  well  as  near  tropical  deep  convection  re- 
10  gions.  The  GW  amplitudes  grow  exponentially  with  height  due  to  the  decrease  in  air 
density.  Different  from  others,  the  GWs  observed  in  AIRS  are  highly  localized,  since 
they  are  dominated  by  high-frequency  components.  This  provides  better  information  on 
GW  sources  from  direct  retrievals  of  GWs,  as  shown  in  Alexander  and  Barnet  (2007). 
The  locations  of  GW  peaks  vary  significantly  from  year  to  year.  The  700  hPa  wind  is 
15  used  to  evaluate  the  GW  generation,  showing  that  stratospheric  GWs  tend  to  grow 
into  great  amplitudes  at  the  places  with  wind  speed  greater  than  10ms_1  collocated 
with  the  upper-level  jet  maxima.  In  the  tropics,  the  deep  convection  controls  the  GW 
variance  peaks. 

The  inferred  preferred  zonal  propagation  directions  of  GWs  seen  in  AIRS  are  oppo- 
20  site  to  the  mean  zonal  wind  directions.  Topographic  GWs  in  the  winter  high-latitudes 
tend  to  propagate  westward  relative  to  the  mean  flow,  while  convectively  excited  GWs  in 
the  tropics  and  subtropics  tend  to  propagate  eastward  relative  to  the  mean  winds.  The 
preferred  wave  propagation  becomes  clearer  with  increasing  altitude,  suggesting  more 
wave  filtering  and  increasing  GWD  on  the  mean  flow.  The  GW  variance  contours  do  not 
25  exactly  follow  the  wind  contours,  suggesting  less  filtering  of  the  high-frequency  GWs. 
Since  the  high-frequency  GWs  tend  to  have  larger  horizontal  phase  speeds,  there  are 
fewer  that  meet  the  critical  levels  in  the  lower  atmosphere  compared  to  lower-frequency 
inertia  GWs.  Therefore,  they  become  more  important  in  the  mesosphere  and  above. 
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An  analogous  approach  of  determining  wave  propagation  direction  has  also  been 
carried  out  on  Aura  MLS  observations  (Wu  and  Eckermann,  2008),  where  wave  pre¬ 
ferred  meridional  propagation  direction  is  inferred  in  a  similar  way  considering  the  fact 
that  MLS  scans  along  the  orbit.  One  can  therefore  obtain  a  more  complete  picture 
5  of  the  GW  horizontal  propagation  direction  from  the  combined  AIRS  and  MLS  results. 
Recent  studies  show  that  convective  GWD  parameterizations  are  sensitive  to  the  wave 
propagation  direction,  which  is  set  to  be  a  free  parameter  in  the  model  (Song  and  Chun, 
2005,  2008).  Our  results  from  AIRS  will  certainly  contribute  to  constraining  GWD  pa¬ 
rameterization  schemes  in  GCMs. 

io  The  new  findings  with  AIRS  GW  variance  differ  from  other  GW  observations  in  three 
aspects: 

Firstly,  AIRS  90-footprint  FOV-depedent  GW  variance  above  mountain  ranges  and 
deep  convection  occurs  with  preferred  angles.  Because  of  the  AIRS  steep  viewing 
angle,  this  finding  indicates  that  really  high-frequency  GWs  with  large  amplitudes  are 
is  present,  and  mainly  excited  from  orographic  and  convective  sources.  It  requires  further 
investigation  on  why  certain  angles  are  preferred,  which  is  a  parameter  missing  in  all 
GWD  parameterization  schemes. 

Secondly,  the  QBO  signals  observed  in  AIRS  are  different  from  other  measurements 
that  mainly  observe  inertia  GWs  and  Kelvin  waves.  With  relatively  fast  vertical  group 
20  velocity,  these  high-frequency  internal  GWs  seen  in  AIRS  are  more  difficult  to  be  filtered 
out  by  QBO  winds,  and  hence  they  are  less  important  in  the  formation  and  descend¬ 
ing  of  QBO  phases.  As  more  of  a  direct  response  to  the  wind  filtering  effect,  AIRS 
GWs  tend  to  propagate  in  opposite  directions  to  QBO  winds.  Since  the  AIRS  GWs  are 
slightly  suppressed  during  the  QBO  westerly  phases,  the  AIRS  GWs  might  be  more 
25  important  to  QBO  westerly  phase  than  to  the  easterly  phase.  The  QBO  has  not  been 
well  simulated  so  far  in  GCMs.  On  one  hand,  the  coarse  vertical  resolution  results 
in  deficiency  in  resolving  some  inertia  GWs.  On  the  other  hand,  the  poorly  observa- 
tionally  constrained  GWD  parameterization  schemes  cannot  precisely  represent  the 
characteristics  of  high-frequency  GWs.  AIRS  not  only  provides  abundant  information 
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on  the  latter  on  to  improve  our  understanding  of  high-frequency  GWs,  but  also  show 
potential  to  observe  low-frequency  inertia  GWs. 

Lastly,  the  most  important  and  encouraging  finding  from  this  work  is  that  AIRS  can 
observe  GWs  with  vertical  wavelength  smaller  than  the  thickness  of  the  WFs.  By  as- 
5  suming  inertia  GWs’  amplitudes  varying  only  with  background  stabilities,  the  variance 
from  the  inertia  GWs  is  measurable  by  AIRS.  AIRS  appears  to  be  highly  selective  on 
high-frequency  GWs,  and  therefore  we  literally  can  separate  the  GW  variance  between 
the  inertia-gravity  and  high-frequency  components.  Moreover,  this  new  capability  will 
help  us  to  further  understand  GW  characteristics  as  observed  by  other  sensors. 

10  This  paper  provides  the  first  climatological  report  of  GW  characteristics  from  nadir 
viewing  sounders.  Since  AIRS  is  sensitive  to  high-frequency  GWs,  the  results  from  this 
study  also  provide  the  first  global  survey  of  the  characteristics  and  properties  of  high- 
frequency  GWs.  Some  large-amplitude  GW  events  as  well  as  FOV-dependent  variance 
can  be  studied  in  further  detail  with  models.  Improving  the  GWD  parameterizations  with 
15  the  AIRS  GW  observations  is  another  important  task  to  explore  in  future  study. 


Appendix  A 

Channel  numbers  and  the  estimated  noise  level 

The  channel  numbers  we  used  in  this  paper  are  listed  below  in  Table  A1.  Together 
20  we  also  give  the  estimated  noise  from  “3-pt  running  average  window”  method  and  the 
provided  instrument  NEdT  on  30  August  2002.  The  minimum  detectable  GW  variances 
for  zonal  means  and  monthly  maps  are  listed  in  the  last  two  columns.  Bold  numbers 
in  the  2nd  column  are  channels  that  locate  at  the  wings  of  the  radiance  spectrum  and 
have  no  paired  channels.  These  channels  are  particularly  sensitive  to  small  drift  in  the 
25  wavelength  while  others  are  stable. 
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Appendix  B 


Some  details  about  the  simulation 


In  the  simulation  mentioned  in  Sect.  4.3,  we  assume  the  imported  wave  amplitude 
varies  with  N.  The  reason  is  as  follows.  Eckermann  (1995)  derived  temperature  vari¬ 
ance  as: 


T'2  =  [i  -  (/y®)2J  (u'2  +  v,2Sj 

Q 


(B1) 


where  g,  u' ,  and  v'  are  the  gravitional  acceleration,  zonal  wind  perturbation,  and  merid¬ 
ional  wind  perturbation  caused  by  GWs,  respectively.  T'2  =  Tl2/T2,  where  T'  is  the  ab¬ 
solute  temperature  perturbation,  while  7"0  is  the  background  temperature.  The  overbar 
denote  the  column  average. 

Think  of  it  physically.  Larger  N 2  corresponds  to  more  stable  atmosphere,  which 
means  the  potential  temperature  increases  more  rapidly  with  height.  Therefore,  same 
oscillations  in  height  would  cause  larger  temperature  perturbations.  By  assuming  that 
waves  neither  break  nor  meet  saturations,  the  vertical  group  velocity  varies  inversely 
with  A/2,  so  that  the  wave  activity  coming  up  from  below  accumulates  in  the  lowermost 
tropical  stratosphere  even  if  the  waves  are  initially  conservative  there. 

The  zonal  mean  values  of  N2  at  January,  2005  are  shown  in  Fig.  1.  The  input 
2-D  wave  amplitude  A  is  assumed  to  be  proportional  to  N  at  the  equator,  which  is 
approximated  as: 


A(x,z) 


-(z-z60) 

40cos(^P;  +  ^P)e  ^  ,  if  z>z60 

Az  Ax  ^  ^ 

A0 cos(^p:  +  ifz<z60 


(B2) 


where  z60  ~  21. 5  km  roughly  corresponds  to  the  geometry  height  at  60hPa.  x  is  the 
along-track  distance,  and  z  is  the  height. 
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The  weighting  function  WF(x, y,z)  for  AIRS  can  be  calculated  directly  from  radiative 
transfer  models  (e.g.,  Eckermann  et  al.,  2007).  Here  we  use  a  2-D  function  combined 
with  two  Gaussian  type  functions  to  approximate  the  actual  weighting  function.  wf(x,z ) 
at  the  nadir  at  80  hPa  is  defined  as: 


5  wf  {x  ,z)  =  • 

(  /-(z-zB0)\2  (  x  \2 

I  e'  V2x6.1  /  *6.75/ 

if  z  >  z, 

|  /-(z-z80)\2  /  x  \2 

^  eV  v^xl  5  /  6'  ^*6.75/ 

if  z  <  z. 

and 

WF(x,z)  = 

wf(x,z ) 

j[wf(x,z)dxdz 

(B3) 


(B4) 


An  example  of  this  weighting  function  has  been  included  in  Fig.  11a.  At  the  outmost 
off-nadir  view,  the  width  of  this  weighting  function  on  along-track  direction  is  broadened 
io  by  a  factor  of  5. 
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Table  A1.  A  list  of  channel  numbers,  noise,  NEdT,  and  minimum  detectable  GW  variance  at 
each  pressure  level.  Bold  numbers  in  the  2nd  column  are  channels  that  are  unpaired.  See 
context  for  details. 


Pressure 

(hPa) 

Channel  numbers 

Noise 

(K2) 

NEdT 

(K2) 

Min.  detectable  GW 
var.  (xIO-3  K2) 
Zonal  Map 

mean 

2 

74 

0.149 

0.165 

3.78 

26.64 

2.5 

75 

0.147 

0.166 

3.72 

26.22 

3 

76 

0.143 

0.161 

3.63 

25.55 

4 

77 

0.145 

0.160 

3.66 

25.80 

7 

78 

0.153 

0.162 

3.88 

27.34 

10 

79 

0.182 

0.172 

4.62 

32.53 

20 

81, 82 

0.084 

0.078 

2.14 

15.05 

30 

102,  108,  114,  120,  125,  126 

0.039 

0.029 

0.98 

6.88 

40 

64,  88,  90,  94,  100,  106,  118 

0.033 

0.028 

0.83 

5.86 

60 

66,  68,  70,  86,  87,  91 , 93,  97,  130 

0.026 

0.018 

0.66 

4.68 

80 

92,  98,  104,  105,  110,  111,  116, 
117,  122,  123,  128,  129,  134,  140 

0.020 

0.011 

0.50 

3.54 

100 

132,  133,  138,  139,  149,  152 

0.026 

0.014 

0.67 

4.73 
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Fig- 1.  Cartoons  showing  the  wave  smearing  effect  from  different  views,  (a)  Westmost  off-nadir 
view,  where  in  this  particular  case  wave  signals  are  largely  smeared  out;  (b)  Eastmost  off-nadir 
view,  where  the  strongest  variance  should  be  observed.  Red  parallel  rectangles  represent  the 
phase  lines  of  the  GW,  while  bold  red  arrows  denote  the  mean  wind  direction. 
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Fig.  2.  AIRS  visibility  as  a  function  of  the  along-track  wavelength  with  3-pt  (black),  7-pt  (blue) 
and  15-pt  (red)  running  smooth  window  applied. 
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(o)  West 


(b)  East 


(c)  West  -  East 


Fig.  3.  Zonal  mean  of  the  GW  variance  at  the  westmost  view  (a),  eastmost  view  (b)  and  the 
difference  between  the  westmost  and  eastmost  views  (c)  as  a  function  of  latitude  and  height 
for  January,  2005.  Monthly  mean  zonal  winds  obtained  from  UK  Met  Office  are  contoured  in 
solid  (westerly)  and  dotted  (easterly)  lines  with  an  interval  of  1 0  m  s_1 .  The  color  scale  is  linear 
by  taking  a  square-root  of  the  values.  In  (c),  positive  (negative)  values  correspond  to  preferred 
eastward  (westward)  propagation  direction. 
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Pressure  (hpa) 


(a)  West  (b)  East 


Fig.  4.  Same  with  Fig.  3,  except  for  July  2005. 
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X  lO-(K')  (b)  East,  2.5  hPa 
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Fig.  5.  Geographical  distribution  of  GW  variance  on  a  2°  x  2°  grid  at  westmost  view  (left  column), 
eastmost  view  (middle  column)  and  the  difference  (right  column)  for  January  2005  at  2.5,  10, 
and  80hPa  from  top  to  bottom.  3-pt  smoothing  has  been  further  applied  to  make  the  signals 
stand  out.  In  the  left  two  columns,  values  smaller  than  2 o  are  uncolored,  where  a  is  the 
minimum  detectable  variance.  In  the  rightmost  column,  absolute  values  smaller  than  V2 a 
are  whitened,  where  \/2a  is  the  smallest  detectable  difference.  The  UK  Met  Office  monthly 
mean  zonal  winds  at  corresponding  levels  are  contoured  in  black  solid  (westward)  and  dotted 
(eastward)  lines.  Total  wind  velocity  greater  than  10  ms-1  areas  at  700  hPa  are  hatched  with 
brown  lines. 
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(a)  West,  2.5  hPa 


X  UT2(K2)  (b)  East,  2.5  hPa 


(h)  West,  80  hPa 


X  10'2(K2)  (j)  East,  80  hPa 


Fig.  6.  Same  as  Fig.  5,  except  for  July  2005. 
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(c)  West  -  East,  2.5  hPa  X 


(f)  West  -  East,  10  hPa  X 


(j)  West  —  East,  80  hPa 
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Fig.  7.  Multi-year  averaged  timeseries  of  GW  variance  for  different  months  at  westmost  view 
(first  column),  eastmost  view  (second  column)  and  the  difference  (third  column)  at  2.5,  10,  40 
and  80  hPa  from  top  to  bottom.  The  UK  Met  Office  zonal  winds  are  overplotted  in  solid  (dashed) 
lines  representing  westerlies  (easterlies)  with  contour  interval  of  10  ms-1. 
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Fig.  8.  (a-c)  Geographical  map  of  GW  variance  at  eastmost  view  at  three  different  locations 
for  January  2005  at  2.5  hPa.  The  GW  variances  within  the  white  box  of  (a-c)  as  a  function  of 
FOV  numbers  are  plotted  in  solid  for  each  location  at  (d-f)  January  2005,  2.5  hPa;  (h-j)  Jan¬ 
uary  2008,  2.5  hPa;  (k-m)  January  2005,  80hPa;  and  (n-p)  January  2008,  80hPa.  Negative 
(positive)  FOV  numbers  correspond  to  west  (east)  views,  and  0  corresponds  to  nadir  view.  The 
linear  fittings  of  each  FOV  curve  are  overplotted  in  dash-dotted  lines. 
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(b)  West  — East  x 
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Fig.  9.  Monthly  mean  timeseries  (annual  cycle  and  linear  trend  removed)  of  the  zonally  av¬ 
eraged  GW  variance  taken  at  the  westmost  view  (a)  and  the  variance  differences  between 
the  westmost  and  the  eastmost  view  (b).  The  vertical  cross-section  is  taken  at  the  equator 
from  80hPa  to  20hPa.  3-month  running  smooth  window  is  applied  to  the  data  to  remove  the 
sub-seasonal  cycle.  UKMO  monthly  mean  zonal  winds  are  overplotted  in  lines  with  contour  in¬ 
tervals  of  5  m  s_1  and  solid  (dotted)  lines  indicating  eastward  (westward)  winds.  Variance  within 
±0.085  x  1 0"3  K2  and  difference  within  ±0.033  x  1 0"3  K2  are  uncolored. 
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Fig.  10.  Zonal  mean  FOV  curves  at  80hPa  for  3-pt  (red)  and  7-pt  (black)  results  at  January 
2005  at  60°  N  (a),  0°  (b)  and  60°  S  (c).  The  horizontal  axes  are  the  number  of  scanning  points 
along  a  FOV  scan,  where  -(+)45  corresponds  to  westmost(eastmost)  view,  while  0  is  at  nadir. 
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Fig.  11.  (a)  Weighting  function  over  nadir  at  80hPa  as  a  function  of  height  and  along-track 
distance  for  the  numerical  simulation.  The  contours  are  0.5,  0.7,  0.9  and  0.99  of  the  peak 
response.  The  outmost  off-nadir  weighting  function  is  4  times  broader  than  that  of  the  nadir, 
(b)  One  example  of  the  original  imported  wave  convolved  with  the  nadir  weighting  function 
(lines)  and  resultant  wave  with  0.1 5n  phase  shift  (shades).  The  wave  has  horizontal  wavelength 
(Ax)  of  100  km  and  vertical  wavelength  (Az)  of  6  km.  Negative  values  are  contoured  with  dash- 
dot  lines  (light  colors),  and  positive  values  are  denoted  with  solid  lines  (dark  colors).  The 
contour  interval  is  0.2  of  the  maximum  and  minimum  amplitudes,  respectively. 
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Fig.  12.  Variance  differences  between  unshifted  wave  and  the  wave  shifted  by  0.1  5ti  for  Xx  = 

700  km  calculation  (a),  and  the  variance  difference  between  nadir  FOV  and  the  outmost  FOV 
(b)  as  functions  of  Az.  The  AIRS  observed  values  (7-pt  smoothing  window  applied)  are  denoted  & 
by  dash  lines.  See  text  for  details. 

03 

o' 

o 

■0 

CD 

"O 

CD 


Gravity  wave 
information  derived 
from  AIRS  radiances 

J.  Gong  et  al. 


Title  Page 


Full  Screen  /  Esc 


Printer-friendly  Version 


Interactive  Discussion 


11737 


ACPD 

11,  11691-11738,  2011 


2.0 
2.5 
3.0 

4.0 

7.0 
10.0 

20.0 

30.0 
40.0 

60.0 

80.0 
100.0 

-90  -60  -30  0  30  60  90 

i  I  i  I  .  1 0"4  (s“2) 

3.7  3.9  4.1  4.3  4.5  4.7  4.9  5.1  5.3  5.5  5.7  5.9  6.0  6.2  6.4 


D 

03 

o 

c 

03 

C/3 

O’ 

o 


T) 


Q) 

"O 

CD 


D 

O)’ 

o 

c 

C/3 

03 

o' 


“0 


Q) 

"O 

CD 


g 

O)' 

O 

c 

03 

03 

O’ 


Tl 

CD 

“O 

CD 


Fig.  B1.  Zonal  mean  of  the  Brunt-Vaisala  frequency  (A/2)  as  a  function  of  latitude  and  height 
for  January  2005,  derived  from  ERA  Interim  dataset. 


g 

c ri 
o 
c 

03 

03 

o’ 


"0 

CD 

"O 

CD 


Gravity  wave 
information  derived 
from  AIRS  radiances 

J.  Gong  et  al. 


Title  Page 


Full  Screen  /  Esc 


Printer-friendly  Version 


Interactive  Discussion 


11738 


